The Asymmetric Thermal Emission of Protoplanetary Disk 
Surrounding HD 142527 Seen by Subaru/COMICS 1 



Hideaki Fujiwara 2 , Mitsuhiko Honda 3 , Hirokazu Kataza 3 , Takuya Yamashita 4,2 , 
Takashi Onaka 2 , Misato Fukagawa 5 ' 6 , Yoshiko K. Okamoto 7 , Takashi Miyata 8 , 
Shigeyuki Sako 8 , Takuya Fujiyoshi 4 , and Itsuki Sakon 2 

S" : 
o ■ 

O ■ ABSTRACT 

<N : 

s : 

Mid-infrared (MIR) images of the Herbig Ae star HD 142527 were obtained at 
\ 18.8 and 24.5 /im with the Subaru/COMICS. Bright extended arc-like emission 

(N \ (outer disk) is recognized at r = (X'85 together with a strong central source 

\ (inner disk) and a gap around r = &!Q in the both images. Thermal emission 

of the eastern side is much brighter than that of the western side in the MIR. 

oo ; 

O \ We estimate the dust size as a few /im from the observed color of the extended 

emission and the distance from the star. The dust temperature T and the optical 

O ! depth t of the MIR emitting dust are also derived from the two images as T = 

82±1 K, t = 0.052±0.001 for the eastern side and T = 85±3 K, r = 0.018±0.001 
for the western side. The observed asymmetry in the brightness can be attributed 
O-i] to the difference in the optical depth of the MIR emitting dust. To account for 

the present observations, we propose an inclined disk model, in which the outer 
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disk is inclined along the east-west direction with the eastern side being in the far 
side and the inner rim of the outer disk in the eatern side is exposed directly to us. 
The proposed model can successfully account for the MIR observations as well 
as near-infrared (NIR) images of the scattering light, in which the asymmetry is 
seen in the opposite sense and the forward scattering light (near side - western 
side) is brighter. 

Subject headings: circumstellar matter — stars: pre-main sequence — planetary 
systems: protoplanetary disks — stars: individual (HD 142527) 



Introduction 



High angular resolution, multiwavelength imaging of protoplanetary disks can provide 
a unique insight into planet formation processes, since detailed structure of a disk should 
bear traces of disk formation history and of young or mature planets. Mid-infrared (MIR) 
observations especially trace the thermal radiation from warm dust, which allow us to ex- 
amine the temperature structure and the dust distribution quantitatively. They also allow 
us to discuss the structure near the central star because the contribution of the photosphere 
is relatively small and coronagraphic observations are generally not required for the MIR. 

HD 142527 (F7IIIe) is classified as a Herbig Ae star as confirmed by its spectrum 
and photometry, which ex hibit the typical characteristics of a pre-main sequence (PMS) 
star (IWaelkens et al.lll996f). Hipparcos measurements suggest its distance of d = 200+4° pc 
( van den Ancker et ali 1998 ). The stellar age and mass are 1 Myr and 2.5 ± 0.3M Q , respec- 
tively, derive d by comparing the po si tion of HD 142527 i n the HR diagram to theoretical 
PMS tracks (Ivan Boekel et al.ll2005l ). iMalfait et al.l (119981 ) reported that infrared excess of 
HD 142527 is comparable to the stellar luminosity, suggesting that a large amount of cold ma- 
terial is present a round the star. Radio observations show the disk mass as Mdi S k = O.15M 
(lAcke et al.ll2004f). which is mor e massive than a typical disk mass by an order of the mag- 
nitude. Ivan Boekel et al. ( 2004 ) reported the central concentration of crystalline silicate in 
the disk. 



Fukagawa et al.l (120061 ) resolved an almost face-on disk around HD 142527 in the near- 



infrared (NIR) scattered light using the Subaru/CIAO. They discovered arc-like components 
facing each other along the east- west direction ( "banana-split structure" ) with one spiral arm 
extending to the north from the western arc in the disk. From this unusual morphology, the 
presence of an unseen eccentric binary and a recent stellar encounter are suggested. However, 
the coronagraphic mask made it difficult to quantitatively discuss the spatial structure in 
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the vicinity of the central star solely from scattered light observations. 

In this letter, we present high spatial resolution MIR images of HD 142527 obtained with 
the Subaru Telescope. The obtained images provide significant information on the nature of 
the protoplanetary disk around HD 142527. 



2. Observations and data analysis 



We observed the He rbig Ae star HP 142527 with the COoled Mid- Infrared Ca mera and 
Spectrometer (COMICS: LKataza et allboodlokamoto et al.lliooil : ISako etal1l2003h mounted 
on the 8.2 m Subaru Telescope on 2004 July 13 and 2005 August 24. Imaging observations 
with the 18.75 /im (AA=0.9 /im) and 24.56 /zm (AA=0.75 jttm) bands were carried out. The 
pixel scale was 0'/13 pixel -1 . To cancel out the background radiation, the secondary mirror 
chopping was used at a frequency of 0.45 Hz for 18.8 /im and 0.46 Hz for 24.5 /im with a 
10" throw. For the Subaru/COMICS, the residual pattern in the chopping subtraction is 
very small and negligible compared to the brightness of the target in most cases. We did 
not apply nodding for the present observation s. We selected stan dard stars (HD 136422 
for 18.8 /xm and HD 146051 for 24.5 /im) from I Cohen et al.l (119991 ) as flux calibrators and 
reference point-spread functions (PSFs). We observed them before or after the observations 
of the target in the same manner. The observation parameters are summarized in Table [TJ 

For the data reduction, we used our own reduction tools and IRAF. The standard 
chopping pair subtraction and the shift-and-add method were employed. We also corrected 
the effect of the different airmass between the object and the stand ard star by estimating the 
difference in atmospheric absorption using the ATRAN software (jLordl Il992l ) . The derived 
total absolute fluxes are 14.6 ± 0.2 Jy and 19.3 ± 0.1 Jy at 18.8 a nd 24.5 /an, respe ctively, 
and they were consistent with the flux of the ISO-SWS spectrum ( ISloan et al.ll2003l ). 



3. Results 

The obtained images are shown in Fig. [T] (left). A prominent extended radiation com- 
ponent is recognized in addition to a central compact component both in the 18.8 and 
24.5 /im images. By comparing with the standard stars, the fluxes of the central component 
of HD 142527 are measured as 10.4 ±0.08 Jy and 6.7 ±0.05 Jy at 18.8 and 24.5 /im, respec- 
tively. The central component is thought to be the thermal radiation from hot dust within 
the inner disk close to the star because the contribution from the photosphere is less than 
0.1 Jy and it is much smaller than the observed emission at these wavelengths. Fig. [2] (top) 
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shows the profile of 24.5 /im image along the east-west direction. While the profile of the 
central component (r < 0'.'5) at 18.8 /im is consistent with the PSF (FWHM= C/52), that at 
24.5 /mi is broader than the PSF (FWHM= 0'.'61). The convolved PSF by the Gaussian of 
FWHM= (X'28 gave the best fit to the profile of the central component at 24.5 /im. This sug- 
gests that the inner disk may be extended and its size can be estimated to be > 80 AU from 
the brightness of the central component at 24.5 jum simply assuming that the disk is a uni 



form b lackbody. Since N-band spectroscopic observations of HD 142527 by Ivan Boekel et al. 



(120041 ) show broad silicate features, /im-sized large grains are dominant in the inner disk. 
From the fluxes at 18.8 and 24.5 /an, the temperature and optical depth of the MIR emitting 
dust within the inner disk are derived as T = 293 K and r = 4.4 x 10 -4 , respectively, by 
assuming a single size distribution of si licate for simplicity and by using the emissivity of 
2.5 /im-sized astronomical silicate dust (ILaor fc Draindll993l ). 



The extended component is thought to be the thermal radiation from the material of 
relatively low temperatures (T < 150 K) in the outer region of the protoplanetary disk since 
it is brighter at 24.5 /iin than at 18.8 /im. In order to examine the extended emission in 
detail, we subtracted the central component. For the 18.8 /im data, we subtracted the scaled 
reference PSF as the central component from the image of HD 142527. For the 24.5 /im data, 
we subtracted the reference PSF that is convolved by a Gaussian of FWHM= 0'.'28 as the 
central component since the central component may be slightly extended (see above). The 
reference PSFs are very stable and the uncertainties of the PSF subtractions are estimated as 
2.6 mJy arcsec -2 at 24.5 /im around the diffraction ring even in the worst case and less than 
the background noise level of 6.7 mJy arcsec -2 at 18.8 /im. The uncertainties are included 
to the errors of the following results. 

After the subtraction of the central component, asymmetric ring-like structures are 
clearly seen (Fig. [T] right). The eastern side is much brighter than the western side both 
at 18.8 and 24.5 /im. We compare the MIR and NIR images in Fig. They were aligned 
with each other by adjusting the stellar position. The stellar position in the MIR image 
was detemined as the peak position of the central component derived by Gaussian fitting 
and its uncertainty is about 0.1 pixel= 0'.'013. The maximum offset between the position of 
the MIR and NIR emission is 0'/013. The observed extended emission arcs are also seen in 



the NIR image as "banana-split structure" (IFukagawa et al.ll2006l ). Gaps in the north and 



south are also clearly seen both in the MIR and NIR images. However, the contrast of the 
two arcs is quite different from the NIR scattered light images, in which the eastern side is 
brighter than the western side. The spiral arm observed in the NIR is not detected in our 
MIR images. 



To discuss the azimuthal variation in the extended emission, we divided the observed 
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extended emission into two regions as the position angle (PA) = — 180° and 180 — 360°; 
eastern and western regions. The radial profiles of the eastern and western sides are shown in 
Fig. [2J The eastern side is about 3 times brighter than the western side. The radial profiles 
are almost Gaussian-like and the peaks are at about r = (X'85 (= 170 AU) for both sides. 
The peak position and profile width (FWHM) in each region derived by Gaussian fitting for 
the 24.5 jum are listed in Table [2j The widths of both arcs are broader than the beam size. 

By integrating the signals within the range of Of! 5 < r < l'!2, we derived the fluxes of 
the eastern and western regions at 18.8 and 24.5 /im and the flux ratio of i 7 i8.8/im/-^24.5/zm 
as a color indicator of the dust emission (Table [2]). We calculated the predicted color profiles 
of the MIR e mitting dust around HP 142527 for the emissivities of various sized astronom- 
ical silicate (ILaor fc Draind Il993l ) and compared them with the observed colors (Fig. H]). 
We assumed L* = 69L m an d T eff = 6250 K for the photospheric spectrum of HD 142527 



(Ivan den Ancker et al.lll998l ). Fig. H] shows that the observed colors of both regions are in 
agreement with the model calculation for 2.5 /zm-sized silicate. Thus we can suggest the 
typical size of the MIR emitting dust within the outer disk of HP 1425 27 as a few /im when 
we assume dust species as amorphous silicate. iFukagawa et al.l (120061 ) also argue that the 
size of grains responsible for the scattering is > 1 /im because the color of scattered light 
from dust around HD 142527 is similar to the color of the star, being in agreement with ours. 
The dust grains around HD 142527 are larger than the grains in the interstellar medium, 
suggesting that the dust grains in the disk have already grown to some extent due to coagula- 
tion of smaller dust particles. We also derived the temperature and optical depth of the MIR 
emitting dust within each region from t he fluxes of 18.8 /im a nd 24.5 /an, assuming a single 
size distribution of 2.5 /xm-sized silicate ( ILaor fc Draindll993l ). We derived the temperatures 
as T = 82 ± 1 K and T = 85 ± 3 K and the optical depths of the MIR emitting dust as 
r = 0.057 ± 0.001 and r = 0.018 ± 0.001 for the eastern and western regions, respectively. 
While the temperatures of the eastern and western regions are in agreement with each other 
within the errors, the optical depth of the eastern region is 3 times larger than that of the 
western region. The asymmetry in the brightness of the extended emission originates from 
the difference of the optical depth between the eastern and western side. 



While the above results are derived on the assumption of 
3OL . 



69L , Ivan Boekel et al. 



( 120051 ) proposed L* ~ 3OL . Taking this luminosity, then we obtain 1.6 /mi as a typical 
size of grains. The effect of the radiation from the inner disk has also been examined by 
assuming that the inner disk emits as a blackbody of T = 1200 K and L ~ 0.5L*, which 
is estimated from the spectral energy distribution (SED). The estimated dust size becomes 
2 /im. Thus we conclude that the presence of /on-sized grains in the disk is secure. 
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4. Discussion 

As a global distribution of the material in the disk around HD 142527, we suggest that 
there are a large amount of cold materials in the outer region of r > 170AU (outer disk), 
while there are quite a few materials present in the inner region. Since we can see a drop of 
radiation at r = C/6 in the images, a gap in the material density must ex ist. The existenc e 



of the material gap is compatible with the MIR flux gap seen in the SED (lAcke et al.l 120041 ). 



Leinert et al.l (120041 ) suggested the presence of "a nearby very red source" whose temperature 
is T ~ 70 K as the origin of the large far-infrared excess seen in the SED. However, our 
observations demonstrate that the large far-infrared excess is attributed to the outer disk 
associated with the star, not to "a nearby very red source" . 

As we mentioned in £|3l the asymmetric brightness of the outer disk of HD 142527 can be 
attributed to the difference in r of the MIR emitting dust. NIR scattered light observations 
of HD 142527 by CIAO show that the western side of the extended disk is brighter than 



the eastern side by 2.0 mag (IFukagawa et al.l 120061 ). 850 /zm continuum observations by 



Submillimeter Array (SMA), which is thought to be optically thin to dust emission and 
a good tracer of the total amount of dust, show an east- west symmetric structure of the 
extended disk (Ohashi et al, in prep.). To explain these observations, we propose a model 
that the dust disk is geometrically thick and inclined as the eastern side is farther to us and 
the inner rim of the eastern side of the outer disk is exposed to us. Fig. [5] shows a schematic 
view of the model of the HD 142527 system. Since the inner edge and rim are hottest within 
the outer disk and affect the MIR brightness most strongly, the emitting region of the outer 
disk appears like a narrow ring. Therefore the optical depth we mentioned above does not 
directly correspond to the vertical depth of the MIR emitting layer of the disk, but rather 
indicates the area size of the MIR emitting region within the beam. If the disk is inclined 
as the eastern side is farther to us, the eastern inner rim of the outer disk can be seen 
directly while the western one cannot because the cooler materials farther out are obscuring 
the inner hot dust. We suggest that the difference of visibility of the inner rim makes the 
strong asymmetric MIR brightness of the outer disk. The model we are pr oposing here can 



also explain that the western side is brighter in the NIR scattered light (IFukagawa et al. 



20061 ). In the scattered light, we observe back scattered light in the eastern side and forward 
scattered light in the western side. For grains those sizes are similar to or larger than the 
wavelength in question, the forward scattering is much stronger than the back scattering, 
resulting in the observed asymmetry in the NIR. An axial symmetric disk inclined to the east- 
west direction is also consistent with the symmetric structure of the extended disk traced by 
850 /mi radiation. It is not necessary to consider an axial asymmetric material distribution 
such as a circumbinary disk to account for the observations. 
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Details of the appearance of a protoplanetary disk, including a "peculiar h orseshoe-like 



featur e" for models with an inclination of 30°, were predicted in simulations of iKessel et al. 



( 119981 ). A model 12 /zm intensity map in their paper shows a morphology strikingly similar 
to our images of the outer disk of HD 142527. NIR interferometric observations of LkHa 101 
show that an asymm etric arc-like radia tion component is present in addition to a central 



compact component (ITuthill et al.ll2002l ). In the case of LkHa 101, the asymmetric arc-like 
structure is thought to be attributed to the hot dust within the inner rim of the inner hole 
made by dust sublimation. Though the observed spatial scale and wavelength are different, 
HD 142527 and LkHa 101 agree on the asymmetric arc-like morphology of thermal emission. 
In the case of HD 142527, photoevaporation of dust can be ruled out as the cause of the gap 
because the position of gap is far from the central star. A secondary object, if any, may be 
able to create the gap. High spatial resolution radio observations will be required in order 
to discuss the origin of the gap in detail. 
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Fig. 1. — Observed images of HD 142527 (left) and central-component-subtracted images of 
HD 142527 (right). These images are in 18.8 /iin (top) and 24.5 jum (bottom). The color 
bars in the right are in the units of Jy arcsec -2 . 
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Distance from the Central Star (arcsec) 



Fig. 2.— 24.5 fim radial profiles of the eastern (PA=60-120°) and western (PA= 240-300°) 
side. Top: pre-subtraction of the central component; Bottom: post-subtraction. The profiles 
of the subtracted component and PSF are also plotted with the dashed and dotted lines, 
respectively. 
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Fig. 3. — Central-component-subtracted 24.5 jum image of HD 142527 (color) and if-band 
scattered light image with CIAO (contour). The color bar in the right is in the units of Jy 
arcsec -2 . 
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Fig. 4. — Color and position of the eastern and western arc. The color profiles calculated by 
using the emissivity of astronomical silicate are also plotted with the solid lines. 
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Fig. 5. — A schematic view of the model of the HD 142527 system. The gap between the 
inner and outer disk exists at r = 80 — 170 AU. 
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Table 1. Summary of the observations of HD 142527 



Object Filter Date[UT] Integ.[sec] AirMass Comments 



HD 142527 Q18.8 2004/7/13 502 2.1 

HD 136422 Q18.8 2004/7/13 201 1.8 standard 

HD 142527 Q24.5 2005/8/24 1404 2.3-2.4 

HD 146051 Q24.5 2005/8/24 201 1.2 standard 
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Table 2. Derived physical parameters for each region 



Region 


F18.8[Jy] 


F24.5[Jy] 


Position 


Width 


F18.8/F24.5 


T[K] 


T24.5 (im 


Center 


10.4 ± 0.08 


6.7 ± 0.05 






1.55 ±0.02 


293 ± 8 


4.4 X 10~ 4 


East 


1.98 ± 0.05 


5.52 ± 0.06 


Of.' 85 


Of.' 87 


0.34 ±0.01 


82 ± 1 


0.057 ±0.001 


West 


0.77 ± 0.05 


2.12 ± 0.06 


0"86 


0'.'94 


0.37 ±0.03 


85 ± 3 


0.018 ± 0.001 



